Distinctness of ensembles having the same density matrix 
and the nature of hquid NMR quantum computing 
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We show that difFerently constructed ensembles having the same density matrix may be physically 
distinguished by observing fluctuations of some observables. An explicit expression for fluctuations 
of an observable in an ensemble is given. This result challenges Peres's fundamental postulate and 
seems to be contrary to the widely-spread belief that ensembles with the same density matrix are 
physically identical. This leads us to suggest that the current liquid NMR quantum computing is 
truly quantum-mechanical in nature. 
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States of an ensemble can be described by a density 
matrix. Unlike a pure state, a mixed state may be pre- 
pared in many different ways. Differently-prepared en- 
sembles may have the same density matrix. It is com- 
monly accepted that ensembles having the same density 
matrix (ESD) can not be distinguished by any conceiv- 
able measurements Q, |M and hence should be con- 
sidered as physically identical. Peres expressed this 
as a postulate: "The p matrix completely specifies all 
the properties of a quantum ensemble." However, some 
other authors expressed their caution. Fano used prepa- 
ration history to specify an ensemble 0- d'Espagnat 
emphatically pointed out that different ensembles hav- 
ing the same density matrix may be physically distin- 
guished, and in Ref. |^, he gave a concrete example for 
a qubit ensemble. Apparently, d'Espagnat's point of view 
has not been widely noticed as many in the community 
tend to believe that ESD can not be physically distin- 
guished. Clarifying this issue has become extremely im- 
portant since in recent years the concept of mixed state 
has been extensively used in quantum computation and 
quantum communication. The prevalent understanding 
on ESD has already led to serious questions about 
the nature of quantum computation: the liquid nuclear 
magnetic resonance quantum computation (NMR QC) is 
not truly quantum-mechanical because of lack of entan- 
glement. 

The question by Braunstein et al. (2| raised upon 
rewriting the effective pure state density matrix in terms 
of a convex decomposition of density matrices of product 
states 



p = (1 - e)Mrf 



(1) 



In Eq. C^l, e is called the polarization and Md is a d x d 
unit matrix divided by d. They then concluded that if 
e in the equation is sufficiently small, the effective pure 
state in NMR QC can be replaced by an ensemble with 
separable states, and therefore, there is no entanglement. 
They thus suggested that the current NMR QC should 



only be considered as a classical simulation but not as a 
true quantum computation. Though Laflamme Q did 
not agree with the conclusion of Braunstein et al. in 
7], he did agree that in each step of the present NMR 
QC experiment the state has no quantum entanglement. 
However, there is one puzzle: if one believes that without 
entanglement, NMR QC is still more powerful than a 
classical computing, then what is the source of the NMR 
QC power, as DiVincenzo questioned 0? 

In Ref. some of us attempted to understand this 
subject. The main points were: 1) The density matrix de- 
scribes the state of "an average particle" in an ensemble. 
It does not describe the state of any individual particle in 
an ensemble; 2) Entanglement is a property of the wave 
function of an individual particle, such as a molecule in 
a liquid NMR sample. Separability of the density matrix 
can not be used to measure entanglement of an ensemble; 
3) The state evolution in a bulk-ensemble NMR compu- 
tation is a quantum process; 4) The coefficient e in Eq. 

is a measure of simultaneity of the molecules in an 
ensemble, which reflects the intensity of the NMR signal 
and has no significance in quantifying the entanglement 
in a bulk ensemble NMR system. Based on these argu- 
ments, it was suggested 10] that the current NMR QC 
is truly quantum-mechanical. 

In this Letter, we extend the study by pointing out 
explicitly that the conclusion that ESD are physically 
identical, which is a form of Peres's postulate and the 
pillar stone of Ref. that has led to doubts on the 
quantum nature of the current NMR computation, is not 
general. We begin with an examination of the mixed 
state concept. Then we focus our discussion on the so- 
called proper mixed state that describes the state of an 
ensemble. We show how ESD can be distinguished by 
considering fluctuations of some observables. We flnally 
come to the question about the nature of quantum com- 
puting, and conclude that the current NMR QC is gen- 
uinely quantum-mechanical. 

A mixed state is the description of state of a quantum 
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system "if we do not even know what state is actually 
present - for example, when several states 4>2t-- with 
respective probabilities wi, W2, ■■■ {wi > 0, W2 > , 
• • •, wi + W2 + • • ■ — 1) constitute the description" 
ll^ . Operationally, a mixed state is usually associated 
with certain averaging processes. There are three ways 
to obtain a mixed state. The first way, which was also 
used by von Neumann in his original papers 11, 12], is 
to average over an ensemble of particles. For simplicity 
we call these particles molecules. Suppose an ensemble 
contains N molecules, with Nwi molecules in state \<j>i), 
Nw2 molecules in state |02), and so on. Then the state 
of an averaged particle in this ensemble is a mixed state, 
described by the density matrix 

p = Wi|0i)(0i| + W2\cj)2){(i)2\ H • 

d'Espagnat referred this kind of mixed state to as a 
proper mixed state 0. 

The second way is to make the averaged state of a 
molecule entangled with other degrees of freedom. By av- 
eraging over the other degrees of freedom using the trace 
operation, the reduced density matrix may become the 
one representing a mixed state. For example, for a quan- 
tum system containing two particles A and B in state 
or equivalently pab = the reduced 

density matrix for particle A is 

PA = TrsipAB)- (2) 

The mixed state arising from this way is called by 
d'Espagnat as an improper mixed state. d'Espagnat 
stressed that an improper mixed state can not be asso- 
ciated with a proper mixed state with the same density 
matrix 6]. 

The third way is to obtain the averaged state of a 
particle over a random variable, or over a period of time. 
For instance, when a spin-i particle is subject to a ran- 
dom kick interaction of the form K = exp[iazO/h], where 
is a random variable with distribution P{0), the av- 
eraged state of this particle over the random variable, 
(P) = /_r P{0)p'd9, normally gives a mixed state, al- 
beit at any instant the particle is in a pure state. 

In the following we consider mainly the proper mixed 
state. The expectation value of observable A for an av- 
eraged molecule in an ensemble is 

(i) = Tr(M), (3) 

and expectation value of A for the whole ensemble of N 
molecules is 

{A)e - N{A) = NTiipA). (4) 

From the expectation values in ||2J) and Q), ESD can of 
course not be distinguished. However, this does not ex- 
clude the possibility of distinguishing them with other 
physical means. d'Espagnat has given an explicit ex- 
ample. Ensemble Sj is prepared by putting N/2 qubits 
in state \z) = |0) and another N/2 qubits in state 



I — z) = |1). In ensemble Su, N/2 particles are put 
in state | + a;) = (|0) + |1))/V2, and another N/2 parti- 
cles in state | — a;) = (|0) — |l))/\/2. The density matrices 
for the two ensembles can be written as 

Si = ^(|0)(0| + |1)(1|), (5) 

Sii - l{\+x){+x\ + \-x){-x\). (6) 

These ensembles are equivalent to those in the example 
used by Peres to show the indistinguishability of ESD 
[T3| . In the matrix form, they both lead to the same 
matrix: M2. However, d'Espagnat used the fluctuation 
of the 3rd component of the spin operator 

N 
i=l 

where the summation is over all molecules in the en- 
semble, to distinguish the two ensembles. For ensemble 
Si, the fluctuation (AS/. = 0, and for ensemble Su, 
(AE//,,)^ = Vn. 

Preskill 14] also noticed that the two ensembles can be 
physically distinguished by using a different method. In 
his example, Alice and Bob share N pairs of particles in 
state IV'As) = ^(|0)yilO)B + ll)yill)B)- Bob measures his 
particles in either cr^ or a^, and Alice's particles is then 
prepared in either Sj or Su. Bob tells Alice of the result 
of his measurement for each particle, but does not tell her 
the measured apparatus. However, Alice can distinguish 
the two ensembles with the information at hand: Alice 
measures each of the particles in her ensemble using the 
(Tz-basis and compares her result with Bob's. If her result 
has a perfect agreement with Bob's, then she knows that 
the ensemble is Sj. Otherwise it is the ensemble Su. 
Hence Alice "does have a way to distinguish Bob's two 
preparation methods" p^ . 

It has become clear from the above discussion that dif- 
ferent preparations of ensembles of qubits can be distin- 
guished physically. We now generalize d'Espagnat's idea 
and show that there exists a general formalism to distin- 
guish ESD. By measuring fluctuations, ensemble prepa- 
ration can be distinguished up to composition: the num- 
ber of particles in each of the pure state in an ensemble. 
We consider different preparations of an ensemble lead- 
ing to the same composition as the same. For instance, 
preparing iVi molecules first in steps of iVi — 3 and then 
adding another 3, or in two equal steps with Ni/2 each, 
will be considered as the same. For an ensemble with N 
molecules with Ni molecules in a state \ipi), the density 
matrix is 

P = ^l^lXV'll + ^IV'2)(^2l + . . . + ^\^n.){i^m\,{8) 

where m is the number of pure states within the ensem- 
ble. Considering an observable 51 for an ensemble 

N 

(9) 

i=l 
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where f2(i) is the observable for molecule i. The expec- 
tation value of for this ensemble can be calculated ac- 
cording to Eq. Q . The fluctuation of is 



m 

\ »=1 



/7VTr(pr!2)-^7V^((V.|l}|7^,))2. (10) 



The term X]j -^^((''/'il^lV'i))^ is sensitive to the ensem- 
ble composition. Here, ^{i) is not restricted to one- 
body operators; for molecules that contain composite 
constituents it can usually be many-body operators. For 
instance, the two-body operator in an ensemble of 
molecules with 7 qubits can be written as 



N 



(11) 



\a<b=l 



By choosing different observables {il,A, •••}, one can 
then distinguish different ESD. 

In principle, the state of an ensemble of N molecules 
with composition ||SJ| can be written as 



1^) = IV-I 



\tpl)Ni |'02>Ari- 



1^: 



2)Ni + N2 



(12) 



though wc usually do not have the detailed information 
about which molecule is in what state. We stress here 
that molecules in an ensemble are classically distinguish- 
able, and that this is different from the state in a Bose- 
Einstein condensate where the particles are indistnguish- 
able. The ensemble results (0J and H10|l can also be cal- 
culated directly by using Eq. H12|l . which leads to the 
same conclusion. 

Now we show that the ensemble of the effective Bell- 
state and its product state decomposition presented in 
Ref. can be physically distinguished. The effective 
Bell-state density matrix is 



Pi = (1 - (-)M4 + f^PBell, 



(13) 



where pseii = (|00) + |11))((00| + (ll|)/2. The effective 
Bell-state describes an ensemble with eN molecules in 
Bell-state (|00) + \U))/V2 and (1 - e)iV/4 molecules in 
each of the calculating-basis states |00), |01), |10) and 
1 11). According to Ref. 7], it can be decomposed into 



P2 



(14) 



where Pi = (1 =p cri)/2 represents a pure state polarized 
or anti-polarized along the the three axes, x, y, and z, re- 
spectively. The two different compositions of the effective 
Bell-state can be distinguished by observing fluctuations 
of the observable 'Szz = Si ^iz(*)'''2z(*)- For the ensem- 
ble (|13|l . the fluctuation is {AT,zz)i,e = &VN, whereas 



for the product state expansion (|14(l . the fluctuation is 

{Al^zz)2.E = ■ Hence the two ensembles are distin- 
guished though they both have the same density matrix. 

The decomposition of the effective Bell-state into a 
product states can not be used to infer the conclusion 
that the effective Bell-state is not entangled because en- 
tanglement is the property of an individual molecule, and 
not of an ensemble. We would like to point out another 
inconsistent problem in Ref. in the description of 
the product state decomposition for NMR QC. In the 
effective pure state NMR QC, the QC computation is 
performed on eN molecules, with each molecule work- 
ing as a single quantum computer. They go through one 
state to another in a computing process. However, in the 
product state decomposition, once an entangling opera- 
tion such as the controUed-NOT gate is performed, the 
number of particles within each composition will change. 

If each molecule in an ensemble is in a mixed state 
itself, even at a given instant, there are fluctuations of 
an observable CI [3i ll^ 



(15) 



In such a case, all ensembles are equivalent and no phys- 
ical means can distinguish them. In fact, there is no 
difference at all among them because the compositions 
of the ensembles are the same, and they all are made up 
by molecules in the same mixed state. This fact has been 
used by some authors to criticize d'Espagnat 17] . To this 
objection, we stress that it is impossible to prepare such 
an ensemble using classical probabilities. For instance, in 
the example given by Peres j the state of a photon is 
prepared according to the result of a coin-tossing, cither 
along z or along —z axis. It is true that each photon 
has 50% probability to be in state |0) or |1). But each 
photon is prepared in a definite state, not in a mixed 
state. Similarly, in the BB84 quantum key distribution 
protocol, the state of a qubit has 25% probability in each 
of the four states \z), \ — z), \x), and | — x). But it is 
certainly in one of the four possible states. This ensem- 
ble has the density matrix M2. The same density matrix 
can be realized by preparing each photon with 50% prob- 
ability in \z) and \ — z), respectively. Though these two 
different preparations result in the same density matrix, 
their security aspect is totally different: the BB84 pro- 
tocol is unconditionally secure whereas the one based on 
the latter is completely insecure. 

One possible way to realize such a scenario is that each 
molecule is actually in an improper state. Suppose that 
we have an ensemble of two-qubits molecules, and the 
two qubits are in an entangled state, say 
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\^Pab)i = ^(|00)as + |11>ab) 



(16) 



If we look at the first qubit in each molecule in the en- 
semble, the resulting state is an ensemble of molecules, 
and each molecule is in an improper mixed state. Then 
in this case, by observing the observables of the qubit A 
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alone one can not distinguish this ensemble from another 
ensemble if the two-qubits molecules are all in the state 

\^AB)2 = ^m)AB + mAB). (17) 

However, we know by now that by looking at the 
fluctuations of a two-body operator, for example, 
'^iLi'^A,z{i)<^B,z{i), one can distinguish between these 
two ensembles. 

Which picture is closer to the liquid NMR QC, an 
ensemble of molecules with each molecule in a definite 
pure state at a given instant, or an ensemble with every 
molecule in the same mixed state? We favor the first 
picture. In the first picture, at a given instant, every 
molecule in the NMR sample is in a pure state and each 
remains in the pure state for some time before changing 
to another pure state. In a long time scale, each molecule 
in an ensemble is in a mixed state. The spin-relaxation 
time Ti in NMR can be viewed as the time period in 
which a molecule remains in a certain quantum state. 
This is supported by the Gorter formula 

where Wm,n is the transition probability rate from level 
m to level n. Hence Ti is the energy weighted time 
in which a molecule remains in a given quantum state. 
A liquid NMR QC is completed within the decoherence 
time, which is much shorter than Ti. During a NMR QC, 
a molecule is well approximated as being in a definite 
quantum state. This picture has been used in the effec- 
tive NMR QC literatures. Gershenfeld and Chuang [T9| 
have explicitly explained this in an effective pure state 
in which there are some numbers of molecules in a def- 
inite quantum state. In addition, nuclear spins are well 



isolated from the environment; its entanglement with the 
environment is weak. For those nuclear spins that are not 
considered in quantum computing, decoupling pulses are 
usually used to disentangle them from the working nu- 
clear spins so that they can not form an entangled state 
that will make the working qubits in an improper mixed 
state. 

In summary, we have shown that for ensembles of 
molecules with the same density matrix, fluctuations can 
be used to distinguish different ensemble compositions. 
This result is a direct generalization of d'Espagnat, and 
it challenges one of Peres' fundamental postulates and 
corrects the common belief that ESD can not be physi- 
cally distinguished. The direct consequence of this result 
is the invalidation of the conclusion in Ref. 0| where 
the authors introduced the decomposition of the effective 
pure state into product states to claim the absence of en- 
tanglement in the current NMR QC experiment. We em- 
phasize that entanglement is the property of individual 
molecules within an ensemble, not of the whole ensemble. 
We suggest that during an NMR QC experiment, quan- 
tum entanglement does exist in individual molecules, and 
this is exactly the source of enhancement of computing 
power in NMR QC. As the only venue capable of demon- 
strating quantum algorithms at present, NMR QC is a 
useful test-bed for quantum computing. Recently, the 
para-hydrogen technique has been used to enhance the 
polarization to near unity [20|. This makes the liquid 
NMR continuing to play the unique role in quantum in- 
formation studies. 
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